Spinodal decomposition of low-density asymmetric nuclear 

matter 

V.Barani'2, M. Colonna\ M. Di Toro^ and A.B. Larionoyi^^ 
^ Laboratorio Nazionale del Sud, 
Via S. Sofia 44, 1-95123 Catania, Italy 
and University of Catania 
^ IFA, Bucharest, Romania and 
^ RRC "I.V. Kurchatov Institute", Moscow 123182, Russia 

We investigate the dynamical properties of asymmetric nuclear matter at low density. The occur- 
rence of new instabilities, that lead the system to a dynamical fragment formation, is illustrated, 
discussing in particular the charge symmetry dependence of the structure of the most important 
unstable modes. We observe that instabilities are reduced by charge asymmetry, leading to larger 
size and time scales in the fragmentation process. Configurations with less asymmetric fragments 
surrounded by a more asymmetric gas are favoured. Interesting variances with respect to a pure 
thermodynamical prediction are revealed, that can be checked experimentally. All these features 
are deeply related to the structure of the symmetry term in the nuclear Equation of State {EOS) 
and could be used to extract information on the low density part of the EOS. 
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I. INTRODUCTION 



During the last years the equation of state {EOS) of nuclear matter (NM) has been studied extensively in the 
symmetric N = Z case. Recent investigations on collisions of radioactive nuclei [11, [3, H, Q and on formation and 
structure of neutron stars [Sj have driven the attention to the properties of strongly asymmetric, N > Z, nuclear 
matter. Hence it appears of relevant interest to investigate equilibrium and non-equilibrium features of asymmetric 
NM and their connection to the used extension of the EOS. 

In this work we study the influence of charge asymmetry on the spinodal decomposition ("5-D") of nuclear matter 
at subsaturation density. The "SD" is the growth of small density perturbations, that leads to liquid-gas phase 
separation, in initially uniform matter located in the low density instability re gion of the EOS phase diagram. In 
symmetric nuclear matter, the " SD^^ was studied extensively in Ref .s ^, 1, 1, HiI, [ll [ll ■ 

To our knowledge, the first publication on the liquid-gas phase separation in asymmetric (neutron star) matter is Ref. 
where the coexistence of asymmetric nuclei with a pure neutron gas was predicted, at zero temperature. Recently, 
the liquid-gas phase transition in asymmetric nuclear matter at finite temperature was studied in Ref.s 
always in an equilibrium thermodynamical approach. For instance, it was demonstrated in Ref. [151]. on the basis of the 
Quantum Statistical Model, that light clusters emitted from neutron-rich systems have, in average, a larger relative 
neutron excess than the initial source. 

Our discussion is more focussed on non-equilibrium properties of asymmetric nuclear matter and, in particular, on 
the possibility to observe collective dynamical formation of clusters, with specific mass and charge contents, on short 
time scales. We expect this mechanism to be important in fast expanding systems, as in collisions of beta-unstable 
nuclei, or as a first step towards preferential equilibrium structures. The latter point, on possible different time scales 
in the clustering process, will be also addressed in this article. 

We consider the unstable growth of density perturbations on the basis of two Vlasov equations, for neutron and 
proton liquids, coupled through the mean field. The influence of the initial asymmetry on the wave length and 
growth time of the unstable modes is studied. In neutron-rich NM the formation of larger fragments is favoured 
and the fragment formation process is delayed in comparison to the case of symmetric NM. The restoration of the 
isotopic symmetry inside the formed heavy fragments, predicted earlier in Refs. [itI . [21|, is clearly observed in 
the " SD" collective dynamical mechanism. Moreover, new non-equilibrium features of the fast clustering processes 
in asymmetric nuclear systems, that cannot be explained within a thermodynamical approach, are discussed. For 
instance, for neutron excess systems the proton fraction in the gas phase is expected to be larger when compared to 
statistical predictions This effect can be related to the "freeze-out" time, preventing the chemical equilibration 
in the dynamics of heavy ion collisions. 

The structure of the article is as follows. In Sect. II a description of the mean- field kinetic approach to the coupled 
neutron and proton liquids is given. We apply the linear response analysis d, H S Elj Il2j for the onset 
of the "SD" instabilities. In parallel, we consider numerical simulations ^2d\ which reproduce all stages of the " SD" 
evolution, including non-linear effects. Results from analitical and numerical solutions of the coupled Vlasov equations 
are presented in Sects. Ill, IV respectively. Summary and conclusions are given in Sect. V. 



II. THEORETICAL APPROACH 



We start from a mean field description of nuclear dynamics based on two Vlasov equations, for neutron and proton 
liquids [l^ [l^, [2lj , coupled through a self-consistent nuclear field : 



dfq{r,p,t) p dfq dUq{r,t)dfq 



dt m dr dr dp 



- . (1) 



Here the subscript q stands for n (neutrons) or p (protons) and /^(r, p,i) is the phase-space distribution function 
(d.f.). For simplicity, in Eq.(l) we neglect effective mass corrections and the difference between neutron and proton 
masses, putting m* — m* — m — 938 MeV . Indeed in the low density region studied here we do not expect to have 
large effective mass corrections. Finally, Uq{r,t) is the self-consistent mean field potential in a Skyrme-like form \2l\ : 

(>Pq \PoJ \PoJ \PoJ 2 dp Po 



where 



2 Po a + 2 2 Po 2 2 
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is the potential energy density; p = + pp and p' = p^ — Pp are respectively the total (isoscalar) and the relative 
(isovector) density; po = 0.16 fm is the nuclear saturation density; Tq = +1 (g = n), -1 {q = p). 

The values of the parameters A = -356.8 MeV, B = 303.9 MeV, a = 1/6 and D = 130 MeV-fm^ are adjusted to 
reproduce the saturation properties of symmetric nuclear matter (binding energy eb — 15.7 MeV/nucleon at p = po, 
zero pressure at p = po, compressibility modulus K — 201 MeV) and the surface energy coefficient in the Weizsacker 
mass formula asurf = 18.6 MeV. We put D' — 40 MeV-fm^ ^ D/3 according to Ref . . that is also close to the 
value D' = 34 MeV-fm^ given by the SKM* interaction [2^. Thus, the term oc D' in the potential energy density ([3]) 
favours the growth of isovector density fluctuations [l3| . 

The potential symmetry energy coefficient is equal to: C(p) = Ci — C2(p/po)", with Ci = 124.9 MeV and C2 = 93.5 
MeV, where the density dependence corresponds to the Skyrme energy density functional of a general kind (c.f. 
Ref.s[2TI. [2^ 23t ). At saturation density the potential symmetry energy coefficient C(po) = 31.4 MeV satisfies the 
condition [24 1 : 
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Cipo) 



where agym — 28 MeV is the symmetry energy coefficient in the Weizsacker mass formula, e^' — 36.9 MeV is the 



Fermi energy for the symmetric system at p = pg. 

As a first step, we apply the linear response analysis to Vlasov Eqs. 
d.f., periodic in time, 6fq{r,p,t) ~ exp(— iwt) we can linearize Eqs. ([T|) : 



For a small amplitude perturbation of the 



iujSfq 



P dSfq 85 fq d6Uq d f^"^ 

m dr dr dp dr dp 



= 



(4) 



where the superscript (0) labels stationary values and SUq is the dynamical component of the mean field potential. 



The unperturbed d.f. /, 



(0) 



is in general a Fermi distribution at finite temperature 



exp (ep - Pq)/T + 1 



(5) 



where — /{2m) + and p,q are respectively energy and chemical potential of the nuclcons of type q. In the 



present work, we neglect finite size effects and consider space-uniform unperturbed d.f. Thus, WrU\ — in Eq. (|4]) 
and we consider plane- wave solutions 5fq oc exp(— iwt -I- ikr). Following a standard Landau procedure 0,|2iIj one can 
derive from Eqs. (|3|) the following system of two equations for neutron and proton density perturbations : 



[1 



+ [K'Xn]Spj, = , 

[l + ^o''Xp]'5pp = 0, 



where 



X<j(w,k) = 



2 dp 



kv 



(0) 



Nq(T) J {2TTh)^ LU- 
IS the long- wave limit of the Lindhard function 0]; v = p/m; 



iO 



kv dcp 



2 dp a/f 



{2TTh)3 



del 



vr 

1-12 



T 



(6) 
(7) 



(8) 



(9) 



is the thermally averaged level density {Nq{0) — mpp^q/ [ii'^fL^), ep^q — pp^^/{2'm), pp^q — h{3n'^ pq)^^^) and, finally 

dU, 



Fr'{k)^Nq,{T) 



5p. 



'11 
'92 



qi = n,_p, q2 = n,_p 



(10) 



are the usual zero-order Landau parameters, where the A:-dependence is caused by the presence of space derivatives 
in the potentials (see Eq.®). For the particular choice of potentials given by Eq.©, the Landau parameters are 
expressed as [2^ : 



Fr'{k)^NqAT) 



Po 



— + {a + l)B^ + Dk^ + (— - D'e)Tq,T, 
' Po 



91 '92 



dC_f[_ 

dp Po " 



— — (^91 +^92) 



d'^C p'- 



dp2 2po 



(11) 
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The dispersion relation connecting eigenfrequencies to wave vectors can be obtained by taking the determinant of the 
system ([5]), © equal to zero : 

(1 + Frx„)(l + F^^^xp) F-^FrXnXp = . (12) 

Since we are interested in the unstable growth, we put in Eq.® lu = ij, where 7 > is the growth rate of the 
instability. Then, in the particular case of T = one can obtain the following simple expression for the Lindhard 
function jS]) (see Ref.[3]) : 

Xqisq) = 1 - Sg arctan(l/s,) , (13) 

where Sq = ^ / [kvp^q)- We stress that without considering the wave number dependence in the Landauparameters 
(fTTjl one has 7 oc fc and, hence, an unphysical growth of the short wave length perturbations is favoured [a. [Til. [T^ . 

Stability conditions of asymmetric nuclear matter against density fluctuations are derived in the Appendix. The 
system becomes unstable if at least one of conditions ([Si]) . ([22]) is violated. 

The linear technique described above is applicable only in the regime of small- amplitude perturbations, and some 
relevant results will be shown in the next section. Then as a second step, see Sect. IV, we consider the more 
powerful numerical solution of the Vlasov Eqs. H]), which is based on the test-particle approach The detailed 

description of the numerical method is given in Refs. [lO, HI]- In this way, we can also have effects from nonlinear 
terms and particle collisions. The latter contribution is indeed always present in the dynamical response of heated 
nuclear matter. However, the main results shown in the next section from a pure mean field approach should not be 
much affected, since a very dilute system is considered. This will be also confirmed from comparisons with numerical 
simulations, where the collision integral is included (see Sect. IV). 



III. UNSTABLE SOLUTIONS OF EXTENDED LANDAU DISPERSION RELATIONS 

We discuss, first, the results of the linear response theory. Fig. 1 shows the instability region (under curves), as 
given by the inequality: 

(1 + i^o"")(l + F^P) - F^^F^'' < , (14) 

with the Landau parameters taken at fc = (see the stability condition ([2^ and text below in the Appendix) in the 
p — T plane for different asymmetries / = {N — Z)/A (a) and in the p — I plane for different temperatures (b). The 
asymmetry leads to shrinking of the spinodal region, reducing both critical temperature and density (Fig. la), in 
agreement with the results of ref. [3| ■ This is indeed a quite general effect due to the attractive neutron-proton effective 
interaction and repulsive neutron- neutron and proton-proton ones 19]: a repulsive symmetry term is softening the 
EOS for asymmetric NM reducing then the low density instability region. 

An increasing temperature also reduces the unstable region in the p — I plane (Fig. lb). 

We have solved the dispersion relation eq. (fT2|) looking at isoscalar growing modes {6pp/6p„ > 0), considering 
various choices of the initial density, temperature and asymmetry of nuclear matter. Fig. 2 reports the growth 
rate F = Im uj{k) as a function of the wave vector k. The growth rate has a maximum Fq = 0.01 ~ 0.03 c/fm 
corresponding to a wave vector value around fco = 0.5 1 fm^^ and becomes equal to zero at fc ~ 1.5fco, due to 
fc-dependence of the Landau parameters, as discussed above. One can see also that instabilities are reduced when 
increasing the temperature. This effect is present also in the symmetric N = Z case jl2l] . 

At larger initial asymmetry the development of the "SD" is slower. One should expect also an increasing of the 
size of the produced fragments due to the " SD" features in asymmetric systems. The effect of the asymmetry on the 
growth time = l/Eo and on the wave length Ao = 27r/fco of the most unstable mode is shown in Fig. 3. It is quite 
clear that the asymmetry dependence of both variables Fq and Aq is more pronounced at higher temperature, when 
the system is closer to the boundary of the spinodal region. 

A better understanding of the "SD" in a two-component system can be achieved by studying the chemical com- 
position of the growing mode. On Fig. 4 we show the asymmetry of the perturbation Ipt — {Spn — 5pp)/{5pn + Spp) 
as a function of the asymmetry of the initially uniform system I = (p^^ — p'^'^)/{pn'^ + Pp^"^)- Without any chemical 
processes, we should expect Ipt — I. However, we obtain Ipt ~ 0.5 /. This means that a growing mode produces more 
symmetric high-density regions (liquid phase) and less symmetric low-density regions (gas phase). Hence, during the 
" SD" , a collective diffusion of protons from low-density regions to high-density regions takes place. 

We see from Fig. 4 that the chemical effect becomes stronger with increasing density. This can be explained by 
the increasing behaviour of the symmetry energy per nucleon with density, in the density region considered here. The 
effect of increasing the temperature goes in the opposite direction, reducing the chemical effect. 
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The conclusion is that the fast "SD " mechanism in a neutron-rich matter will dynamically form more symmetric 
fragments surrounded by a less symmetric gas. Some recent experimental observations from fragmentation reactions 
with neutron rich nuclei seem to be in agreement with this result on the fragment isotopic content : ne arly symmetric 
Intermediate Mass Fragments (IMF) have been detected surrounded by very neutron- rich light ions [2^ . 



IV. NUMERICAL RESULTS: HEATED NUCLEAR MATTER IN A BOX 



The previous analytical study is restricted to the onset of the "S*!?", in a linearized approach. Then numerical 
calculations were performed in order to study all stages of the fragment formation process. In the numerical approach 
we consider nuclear matter in a cubic box of size L imposing periodic boundary conditions. We follow a phase-space 
test particle method to solve the Landau- Vlasov dynamics, using gaussian wave packets [H, [l^, H^. The dynamics 
of nucleon-nucleon collisions is included by solving the Boltzmann-Nordheim collision integral using a Monte-Carlo 
method [27l |. 

We choose the width of the gaussians in order to correctly reproduce the surface energy value in finite systems. In 
this way a cut-off appears in the short wavelength unstable modes, preventing the formation of too small, unphysical, 
clusters [l^l ■ In order to have a correct mean field treatment also on the edges of the box, we have used a self-consistent 
stuffing method, filling a layer of 6fm thickness all around the box with test particles having symmetric positions 
with respect to those on the opposite side in the box. 

The calculations were performed using 80 gaussians per nucleon and the number of nucleons inside the box was 
fixed in order to reach the initial uniform density value. An initial temperature is introduced by distributing the test 
particle momenta according to Fermi d.f. (Eq. (5)). We have checked that we reproduce at equilibrium the right 
"EOS" corresponding to the used effective forces (see Sect. II). 

We have followed the space-time evolution of test-particles in the box with side L — 24 fm for three values of the 
initial asymmetry / = 0, 0.25 and 0.5, at initial density — OApo and temperature T = 5 MeV. The initial density 
perturbation was created automatically due to the random choice of test-particle positions. Results for the initial 
asymmetries / = and / = 0.5, are reported in Fig. 5, (a) and (b) respectively. Figure 5 shows density distributions 
in the plane z = 0, which contains the center of the box, at three time steps t — 0, 100 and 200 fm/c, corresponding 
respectively to initial conditions, intermediate and final stages of the " SD" . Clearly, the growth of the small initial 
density perturbations takes place. 

We have compared the dynamical evolution, as given by the test particle method, with the analytical predictions 
of Sect. Ill (Figs. 2-4). To do this, two variables were constructed: the total density variance (see [9|]) 

a=< (p-p("))2 (15) 

and the correlation function between proton and neutron density perturbations, normalized to the neutron density 
variance, 

p < jPp- Pp^)iPn- Pn'') >all 

-"^P" : rmTT • 'y^'^) 



<{Pn-pl?r >n 

In Eqs. ([T5|) ■ (fT6|) < ... >au denotes the average over all test particles, while < ... >„ denotes the average over neutrons 
only. The densities p, p„ and pp were calculated in the position of the test particle considered by taking contributions 
from gaussians of all test particles. For a dominant plane-wave perturbation we have: 

a oc exp(2rt). 



Opn 



Fig. 6 shows the evolution of a (a) and of the (test-particle) perturbation asymmetry Ip^ — (I — -Rp„)/(1 -I- i?p„) (b), 
for the same initial conditions discussed above, i.e. T = 5 MeV, = OApo and asymmetries / — 0.0, 0.25, 0.5. 

A general feature is the clear linear increase of ln{a) in the time interval 50 < t < 150 fm/c. During the first 
50 fm/c the system is quickly "self-organizing" selecting the most unstable normal mode. Afterwards the variance 
(Eq. HI])) increases exponentially with a time scale given by F = Im uj{k). In correspondence (see Fig. 6b), the 
perturbation asymmetry Ipt reveals also a quick saturation at i ~ 50 fm/c. At earlier times the proton and neutron 
density perturbations are not correlated, but at i > 50 fm/c the correlation of plane-wave type {Spp/Spn = const > 0) 
develops. 

We notice that the time scales necessary to reach the asymmetry value characteristic of the most important growing 
modes, that are quite short in our calculations, generally depend on the structure of the initial noise put in the neutron 
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and proton densities. In our calculations all modes are nearly equally excited. This causes the quick appearance of the 
features associated with the dominant mode. In agreement with analytical calculations, the instability grows slower 
in the case of larger asymmetry. 

For an initial asymmetry I = 0.5, the extracted values of growth time T ~ 0.01 c/fm and perturbation asymmetry 
Ipt ~ 0.24 (see Fig. 6), and of wave length A ~ 12 fm (from the distance between the density distribution maxima in 
Fig. 5b), are in good agreement with the analytical results presented in Figs. 3,4. 

The "S"!)" leads to a fast formation of the liquid (high density) and gaseous (low density) phases in the matter. 
Indeed this dynamical mechanism of clustering will roughly end when the variance (Eq. (jlSp ) saturates [30], i.e. 
around 250 fm/c in the asymmetric cases (see Fig. 6a). 

We will discuss in the following the "chemistry" of the liquid phase formation. In Fig. 7 we report the time 
evolution of neutron (thick histogram in Fig. 7a) and proton (thin histogram in Fig. 7a) abundances and asymmetry 
(Fig. 7b) in various density bins. The dashed line respectively shows the initial uniform density value p ~ 0.4po (Fig- 
7a) and the initial asymmetry / = 0.5 (Fig. 7b). The drive to higher density regions is clearly different for neutrons 
and protons: at the end of the dynamical clustering mechanism we have very different asymmetries in the liquid and 
gas phases (see the panel at 250/ m/c in Fig. 7b). 



It was shown in Refs. [ij, on the basis of thermodynamics, that the two phases should have different 

asymmetries, namely, Igas > Iiiquid, and actually a pure neutron gas was predicted at zero temperature if the initial 
global asymmetry is large enough (J > 0.4) [l4| . Here we are studying this chemical effect in a non-equilibrium 
clustering process, on very short time scales. The interest is in the observation of new features, not expected in a 
thermodynamical picture. 

In our numerical model, we divide the system into liquid-like and gas-like phases as follows. The i-th test particle 
belongs to the hquid- (gas-) like phase if pi > (<) p{t = 0), where pi is the total density in the position of the i-th 
test particle, p{t = 0) is the initial density [sij . 

Fig. 8 presents the time evolution of the asymmetry in liquid and in gas. Strong fluctuations take place at the 
beginning of the evolution. However, at t > 50 fm/c, results are quite stable. In the symmetric case (/ = 0), both 
liquid and gas keep the symmetry in the course of the time evolution. However for non-zero initial asymmetry liquid 
becomes more symmetric and gas less symmetric (neutron-rich) as time goes. The asymmetry in liquid and gas phases 
saturates rather early, at t = 200 ^ 300 fm/c, just on the time scale of dynamical clusterization. 

This suggests an evidence of a two-stage fragment formation process. On the first stage, the fast " SD" takes place 
(t < 250 fm/c, see Fig. 6a). At the end of this stage, due to reduced asymmetry of the perturbation, see Fig. 4, liquid 
acquires a lower asymmetry and gas acquires a higher one. On the second stage {t > 250 fm/c) we have a statistical 
nucleation process: protons and neutrons diffuse very slowly from gas to liquid, see Fig. 9a, with approximately 
equal rates (Fig. 9b). Since the gas is already highly asymmetric and the liquid is more close to symmetry, this 
diffusion process leads to a further growth of asymmetry in the gas and does not change the asymmetry of the liquid. 
Eventually, on long time scales, we can expect to reach the thermodynamical limit of Refs. i-^- almost 

pure neutron gas. 

From our dynamical analysis we can conclude that the second slow mechanism for liquid formation is certainly 
present in a confined system and will be stopped in an expanding case. We deduce that the isotopic contents of the 
gas phase could give a measure of the "freeze-out" time: a very neutron-rich gas will correspond to a slowly expanding 
nuclear system. 

This effect is quantitatively shown in Fig. 10 where we present asymmetries of the two phases at different time stages, 
as given by the test-particle simulations. We choose two " freeze-ouf times, one corresponding to the intermediate 
stage of the " SD" {t = 150 fm/c) and the second inside the slow diffusion process {t = 400 fm/c). The first choice of 
freeze-out time is close to the time scale of the fast multifragment breakup in intermediate energy heavy ion collisions 



The asymmetry of the liquid is, practically, not dependend on the choice of " freeze- out" time. On the contrary, the 
asymmetry of the gas and the number of gas particles are very sensitive to this choice. 



An investigation of fragment formation through Spinodal Decomposition in low-density asymmetric nuclear matter 
has been performed. Important information on the early evolution of the unstable modes is obtained performing an 
analytical linear analysis to the Vlasov equations, for neutron and proton liquids, coupled through the mean field. 
Then a numerical study of nuclear systems with periodic boundary conditions has been applied to describe all stages 
of the fragment formation process. 

We have shown that charge asymmetry changes time and space scales of the fragment formation process. This 
effect could be observed experimentally. 




(c.f. [Tl[l,il). 



V. 



SUMMARY AND CONCLUSIONS 
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With respect to the fragment isotopic distribution, we predict that the " SD" in asymmetric systems is accompanied 
by a coUective diffusion of protons from low to high density regions. Thus the produced heavy fragments will be more 
symmetric than the initial uniform system. This qualitatively corroborates earlier thermodynamical studies [13] ■ 
However the chemical composition of fragments is established just after the finishing of the fast " SD" stage. Hence 
it is essentially defined by the non-equilibrium mean field dynamics and a thermodynamical approach gives only a 
rough schematic description of this process, with a slow nucleation component likely absent in fragmentation reactions 
with radioactive beams. In general in a dynamical multifragmentation process with neutron rich systems we expect 
a proton fraction in the gas phase, i.e. among emitted nucleons and light clusters, larger than the thermodynamical 
prediction. 

Finally we think that different choices of effective interaction, (eg. with different density dependence in the symmetry 
energy coefficient C{p)) will support the qualitative conclusions of our work, since all realistic effective interactions give 
a similar behaviour for the potential symmetry energy per nucleon in the region of subnuclear densities (ssj . However 
the size of the instability region and the isotopic structure of the most unstable collective modes will be certainly 
dependent on the used effective force. Therefore, considering fragmentation reactions with radioactive beams, it would 
be possible to extract important information on the low density part of the EOS in asymmetric nuclear matter, of 
large interest for the understanding of life and structure of neutron stars. 
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Appendix: Stability conditions 

Here we derive the stability conditions of asymmetric nuclear matter against density perturbations. Let us write 
down the variation of the free energy density J- = e — Ta, where e and a are respectively energy and entropy densities, 
over proton and neutron densities, up to the second order [l4| keeping constant temperature and volume : 



ST = UnSpn + fJ'pSpp + 



d^in Spl dup Spl 



dpnOpp 

OPn 



(17) 



dpn 2 dpp 2 

First two terms disappear after integration over the volume, since the total neutron and proton numbers are conserved : 



d^rST : 



dpn 



5pl 



dpp 



Spl 



2^^6pnSpp 

OPn 



(18) 



In order the system to be stable against density fluctuations, the quadratic form in square brackets of eq. (jlSp must 
be positive defined. Thus, we obtain two stability conditions ; 



^^">0 or 5^>0 



dpn dpp 
dp.n dp,p _ /djlp^ 2 

dpn dpp \dpn 



> . 



(19) 
(20) 



The system is unstable if at least one condition is violated. We notice that it is enough to satisfy only one inequality 
in (19), since the other one will be satisfied automatically, if condition (20) is fulfilled. In an equivalent way, the 
stability conditions and (|^D|) can be expressed in terms of the Landau parameters as follows : 



H 
(1 



Fq" " > 



or 1 



F^'PFP"" > 



(21) 
(22) 



where we have used the relation 



dPq' 



^0 ' 



q = n,p 



(23) 
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with Sqqi = 1 a q = q' and 5qqi = Q ii q ^ q' . The relation (j23|) is obtained by considering the variation over 5pqi of 
both sides of the equation: 



P.= 7^.^(4)- (24) 



In the case of the interaction used here (see Eq. ([TT]) ). the condition (PT|) is always satisfied. 

The stability condition pO)) was derived first in Ref. [M]- Then it was shown in Ref. that this stability 

condition can be rewritten equivalently as: 



a I \ I >0, (25) 

"P J T,y \"y J T,P 

where y = Pp/ p is the proton fraction, P = p — J- is the pressure. 

The boundary of the instability region can be obtained just putting = in the Lindhard function ([8]). Since 
Xq{oj = 0, k) = 1, we have from dispersion relation : 

(1 + Fo"") (1 + PPP) ~ F^^F^-^ = . (26) 

In the case of symmetric nuclear matter {F^^ = F^^ , Fq^ — Fq") the stability conditions (|2ip . (P^ are equivalent to 
the system of Pomeranchuk stability criteria Q for isoscalar and isovector motions : 

1 + Fo > , (27) 
1 + J^o > , (28) 

where Fq = F^^ + F^^ and Fg = FJ*" — F^^ are the isoscalar and the isovector Landau parameters [13] ■ 
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Figure captions 

Fig. 1 Spinodal boundaries in density-temperature plane at different asymmetries (a) and in density-asymmetry 
plane at different temperatures (b). Instability regions are under curves. 

Fig. 2 Growth rate of instability as a function of the wave vector, as calculated from the dispersion relation p2|l for 
p = 0.2po, r = 5 MeV (short-dashed hues); p = 0.4po, T = 5 MeV (soHd lines); and p = 0.4po, T = 10 MeV 
(long-dashed lines). Lines are labeled with the asymmetry value /. 

Fig. 3 Wave length Aq (a) and growth time (b) of the most unstable mode as a function of the asymmetry / for 
p^"^ — 0.2po (dashed lines) and p*^"^ = 0.4po (solid lines). Lines are labeled with the temperature T in MeV. 

Fig. 4 Perturbation asymmetry Ipt versus initial asymmetry / for different initial densities and temperatures : p^"-' — 
OApo and T = 5 MeV (sohd line), p^°'> = 0.2po and T = 5 MeV (short-dashed line), p^°'> = 0.2po and T = 10 
MeV (long-dashed line). 

Fig. 5 Time evolution of the density p{x, y) in the plane z = as given by the test particle code for initial density 
— 0.4/30, at temperature T = 5 MeV and asymmetries / = (a) and / — 0.5 (b). Upper panels show 
contour plots of the function p{x, y) and lower panels report the corresponding two-dimensional surfaces. The 
density is in units of fm^^ x lO'^. 

Fig. 6 Time dependence of the density variance (a) and of perturbation asymmetry (b) (see text for definitions) in 
test particle simulations at initial density p*^°) — 0.4po and temperature T = 5 MeV for initial asymmetries 
1 = 0, 0.25 and 0.5 (solid, long- and short-dashed lines respectively). The straight lines on Fig. 6a show linear 
fits to the initial stage of the SD. 

Fig. 7 Time evolution of neutron (thick solid lines) and proton (thin solid lines) abundances (a) and asymmetry (b) 
in different density bins. Calculations refer to the case with initial temperature T = 5 MeV. Initial values of 
density p'"' = 0.4po and asymmetry / = 0.5 are indicated by dashed lines. 

Fig. 8 Time evolution of asymmetries in liquid (solid lines) and in gas (long-dashed lines) for initial asymmetries 
(shown by short-dashed lines) / = 0, 0.25 (a) and / = 0, 0.5 (b). 

Fig. 9 Ratio of total number of particles in gas to total number of particles in liquid (a) and number of neutrons and 
protons in gas (b) versus time. 

Fig. 10 Asymmetries in gas (squares) and in liquid (circles) as functions of the initial asymmetry. Long-dashed 
(solid) lines correspond to the " freeze-ouf time 400 (150) fm/c. In the liquid phase there is no appreciable 
difference between the two lines. 



8X0 0X0 BOO 900 tOO 200 00 




92' T 



OOT 9^0 090 9Z0 




000 
000 



— TOO 



Am 01 =L 

odro=cf 



Am 9=1 ^dro=d 



Am 9=1 ^dz'0=d 



J I I I I I I I I I I I I I I I I I I I I I I I I ^O'O 



— 200 



— eoo 



90 to eo Z'O 



TO 



"1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




00 
00 



— TO 



— Z'O 



Am 01=1 ^dz'o=d 

km 9=1 ^dz'0 = d 



— eo 



to 



1=; 

<rt- 
I— 

O 
P 



CD 
<rt- 



90 




100 200 300 400 500 100 200 300 400 500 

time (fm/c) time (fm/c) 



t=0 fm/c 



t=50 fm/c 



t=100 fm/c t=150 fm/c 



t=250 



150 



100 - 



ji 50 - 




150 



100 - 




50 - 



P (fm-^) 



1.0 



0.8 

0.6 E- 

0.4 r 

0.2 r 

0.0 
0.0 



L(b) 



1.0 
0.8 F- 
0.6 r 



0.4 r 
0.2 r 



0.1 



0.0 



0.0 



0.1 




